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Two Fokker actions and corresponding equations of motion are obtained for two point particles 
in a post-Minkowski framework, in which the field of each particle is given by the half-retarded 
+ half-advanced solution to the linearized Einstein equations. The first action is parametrization 
invariant, the second a generalization of the afhnely parametrized quadratic action for a relativistic 
particle. Expressions for a conserved 4-momentum and angular momentum tensor are obtained 
in terms of the particles' trajectories in this post-Minkowski approximation. A formal solution to 
the equations of motion is found for a binary system with circular orbits. For a bound system of 
this kind, the post-Minkowski solution is a toy model that omits nonlinear terms of relevant post- 
Newtonian order; and we also obtain a Fokker action that is accurate to first post-Newtonian order, 
by adding to the post-Minkowski action a term cubic in the particle masses. Curiously, the conserved 
energy and angular momentum associated with the Fokker action are each finite and well-defined 
for this bound 2-particle system despite the fact that the total energy and angular momentum of the 
radiation field diverge. Corresponding solutions and conserved quantities are found for two scalar 
charges (for electromagnetic charges we exhibit the solution found by Schild). For a broad class of 
parametrization-invariant Fokker actions and for the afnnely parametrized action, binary systems 
with circular orbits satisfy the relation dE = ildL (a form of the first law of thermodynamics), 
relating the energy, angular velocity and angular momentum of nearby equilibrium configurations. 



I. INTRODUCTION 

In constructing numerical initial data for compact binary systems, it may be useful to find exact solutions to the 
Einstein equations with helical symmetry, solutions that are stationary in a corotating frame 

In the case of two charged particles, an exact solution of this kind was found by Schild[ll(: The force on each particle 
is associated with the half-retarded plus half-advanced field of the other. (An earlier but less explicit version is given 
by Schonberg|l2j.) We obtain here an analogous solution for point masses in a post-Minkowski framework. 

A time-symmetric interaction allows one to derive the equations of motion for point-particles from a single action 
integral that is written solely in terms of the dynamical variables of each particle (without mediating field variables) . 
This kind of action at a distance theory has been formulated by Fokker and by Wheeler and Feynman; for their 
treatments of the electromagnetic field, see 0|. We exhibit an analogous Fokker action for linearized gravity. 
The equations of motion are those of a post-Minkowski approximation, in which, as in the electromagnetic case, each 
particle moves in the half- advanced plus half- retarded field of the other. 

A Fokker action is not a true action: Its variation leads to equations in which the endpoints of the action integral 
explicitly appear in integrals that should yield the field of each particle. Only by taking a limit of the varied action as 
the endpoints go to infinity does one recover the correct equations of motion. Invariance of the Fokker action under 
Poincare transformations nevertheless leads to expressions for a conserved energy and angular momentum of an n- 
particle system. For helically symmetric binary systems (in, e.g., models with scalar, electromagnetic, or gravitational 
interactions), the energy and angular momentum of the field is infinite, because the particles have radiated for infinite 
time in past and future. The conserved energy-momentum and angular momentum associated with the Fokker action, 
however, are finite. For circular orbits of particles of mass m and to, we find the surprising result that the form of 
this energy, 

£=^ + ?, (1) 
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is identical for scalar, electromagnetic, and linearized gravitational interactions, described by a parametrization in- 
variant Fokker action. The angular momentum is in each case proportional to the value of the interaction field dotted 
on each free index into the helical Killing vector. We show that nearby circular orbits satisfy the familiar first law, 
dE = ndL. 

Although one can construct formal solutions to the post-Minkowski equations in which the source is gravitationally 
bound, care is needed in their interpretation. A conceptual problem is related to the fact that, at zeroth order in the 
perturbed metric, particles move on geodesies of flat space, not in bound orbits. Because the matter density vanishes 
at zeroth order, the stress-energy tensor T Q/3 = pu a u@ has at first order the form ST a P = (Sp)u a u^ , where u a is the 
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zeroth-order velocity field describing straight-line motion. If one then finds the metric from the perturbed equation 

= -25 [G Q/3 - %TrT af3 } = D(h a/3 - \rf !i h) + 16n(6p)u a u p , (2) 

the perturbed metric will not yield perturbed bound trajectories. 

A resolution to this problem is to look at families of solutions, g a p(s),T a P(s) that describe finite fluid masses in 
bound orbits for all nonzero values of s, despite the fact that s — is flat space. An example can be easily understood 
by first looking at Newtonian gravity, with a family of solutions corresponding to the scaling v ^ e, m/r ^ e 2 of the 
post-Newtonian approximation: That is, consider a family of solutions for two point masses m — m^'s and fh — fh^s 
in circular orbits of radii a and a (independent of s), with speeds given approximately by v^y/s and v^yfs. At 
s = 0, the solution is flat space with no matter. For each nonzero s the solution describes masses in circular orbit, 
and as s — > 0, the period of the orbit increases without bound. An exact solution is given for each mass by equations 
in which the perturbed particle trajectory is found self-consistently. That is, with $ = s$^\ the center-of-mass 
trajectories r and r satisfy 

V 2 $ (1) (r) = 47rm (1) (5(r- r(t)), f = -V$, 

V 2 $ (1) (r) = 47rm (1) (5(f - r(i)), f = -V$. (3) 

Here $(r) is the potential at r due to the mass m. Because v — > as s — > 0, the orbit remains close to a straight line 
for increasingly long times, but for each finite s the orbit is circular . 

In general relativity, one can presumably construct a similar family of solutions with bound orbits whose energy- 
momentum tensor, T a @(s) — e(s)u a (s)u@ (s) + p(s)[g a ^ (s) + u a (s)u^(s)], and metric g a j3{s) are pointwise continuous 
in the parameter s and for which one has flat, empty space at s = 0, 

g a p(Q)=Vap, e(0)=p(0) = 0. 
The equations for the first-order metric in a radiation gauge, 

- 2G« = U{h aP - \^)h = -167tfg, (4) 

have as their source a stress-energy tensor constructed from the perturbed velocity field and density and from the 
unperturbed metric. With the masses shrinking to zero as s — > 0, the motion of each mass is given to linear order by 
the linear field of the other: The self-force at linear order serves only to renormalize the mass. 

The resulting "post-Minkowski" solution in which two masses move in bound orbit, each responding to the linear 
field of the other has the following features: 

(i) It is correct to Newtonian order. 

(ii) The radiation field of the linearized metric is correct to lowest nonvanishing post-Newtonian order (2 1/2 post- 
Newtonian order). 

(iii) The orbit is not correct to first post-Newtonian order. 

Using the linearized metric leads to equations of motion that are missing a term quadratic in the particle masses that 
enters the first post-Newtonian equations. In this sense, the first post-Minkowski approximation for bound orbits 
is a toy model, keeping terms of all orders in v but discarding nonlinear terms in h a p that, for bound states, give 
corrections of order v 2 to the orbit. 

Following recent codes that obtain helically symmetric solutions to nonlinear wave equations 0, 0] , Uryu 0, 0] 
has obtained a neutron star code that solves the full Einstein-perfect fluid equations on an asymptotically flat initial 
hypersurface S, obtaining a solution with exact helical symmetry in the near zone. By solving the full Einstein equation, 
one expects more accurately to enforce the balance of gravity and centripetal acceleration in circular motion. The 
error associated with ignoring the radial motion associated with radiation reaction, however, remains. In future work, 
we anticipate using the point-particle model developed here to estimate the accuracy of such helically symmetric 
initial data sets and corresponding quasiequilibrium sequences, by comparing outgoing point-particle solution (in a 
post-Minkowski framework) to helically symmetric models and sequences. A further problem is related to the fact 
that the codes mentioned above typically diverge when helical symmetry is enforced in a region larger than the near 
zone (beyond one wavelength). Extending the present work to second order in the post-Minkowski approximation may 
be useful in developing codes with larger domains of convergence. Finally, existence of helically symmetric binaries 
is unproved and the second-order (nonlinear) extension may also help in understanding existence and asymptotic 
behavior of models of helically symmetric binaries. 

In Sect.^ we review Fokker actions, obtaining the equations of motion and the form of the conserved 4- momentum 
and corresponding angular momentum. In Sect. IIII Al we introduce a parametrization-invariant Fokker action that 
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describes point particles in the first post-Minkowski approximation. The formalism of the previous section is used 
to obtain equations of motion and conserved quantities; derivations of the explicit forms of the conserved quantities 

are relegated to Appendix ID 31 Sect. IIIIBl thcn introduces the Fokker analog of the quadratic action, / dr—x a x a , 

for affinely parametrized particles, again obtaining equations of motion, momentum, and angular momentum. In 
Sect. ITTTO and Till Dl we consider two particles of masses m and to in circular orbit, finding the equations governing 
the orbit, and computing the system's energy and angular momentum. The corresponding conserved quantities for 
scalar and electromagnetic interactions are summarized. Sect. Hvl is devoted to proving the relation dE = QdL for a 
class of parametrization- invariant Fokker actions and the affinely parametrized action. In Sect. we introduce two 
forms of a post-Newtonian correction term to make the post-Minkowski action accurate to first post-Newtonian order. 
We obtain the corresponding corrections to the conserved momentum and angular momentum. Sect. IVII displays the 
results of a numerical solution to the post-Minkowski equations of motion, with and without the added post-Newtonian 
term. Sect. IVIT1 discusses features of conserved quantities associated with Fokker actions, proposing, in particular, an 
explanation for their finite behavior when the field energy is infinite. A description of interacting scalar charges and 
many of the details of our calculations are presented in appendices. 



II. ACTION AT A DISTANCE THEORY 



An action-at-a-distance theory of interacting classical point charges was formulated by Fokker and by Wheeler and 
Fevnman.|l3l IT^|. In these treatments, one obtains the equations of motion by varying an action integral that is a 
function only of the trajectories (worldlines) of the charges. The price one pays for eliminating the electromagnetic 
field is that the action integral is not a genuine action: One must vary an integral / that involves only a finite segment 
of each particle's trajectory, and the equations of motion emerge from the limit of SI as the trajectories are extended 
to infinity. In other words, the limit must be taken after the variation of the action integral. In addition, the action 
is not an integral over a single parameter time, but instead involves integrals over parameter times associated with 
each particle. 1 We will call an action integral having these properties a Fokker action. 

In Dettman and Schild [l8|, a derivation of the equation of motion, as well as expressions for the conserved 
quantities, the energy momentum and the angular momentum, is presented for a generic Fokker action that may 
include self-action terms and variable mass parameters. 

In this paper, we obtain a similar action for two self-gravitating point particles, in a post-Minkowski approximation. 
We describe the particles by their constant masses to, to and by their trajectories written in terms of position vectors 
a; a (r), x a (f) on a flat background spacetime, with arbitrary parameter times r and f. In this approximation, the 
motion of each particle is determined by a gravitational interaction with the other particle, and there is no self- 
interaction term (any contribution from gravitational self-energy is accounted for in the mass of each particle). The 
two-particle system is then described by a Fokker action of the form 

I{ T i, T 2, Ti, T2) = — to / dr(— x a x a ) 3 — to / df(—x a x a )^ + / dr df A(x — x, x, x), (5) 

J Ti J fi J 7"i J T\ 

where the Fokker analog A of an interaction Lagrangian has the property, 

A(x — x, x, x) = A(x — x, x, x), (6) 

and where 

x a :=^, and x a := ^f-. (7) 
dr dr 

For given values of r, f, A(r, f) is a scalar constructed only from r\ a p,, x a (r), i Q (r), x a (f ), x a {f). It follows that A and 
/ are Poincare invariant - invariant under simultaneous Poincare transformations of the particles' paths. 

As noted earlier, to obtain the equations of motion, one must compute the variation of the action integral before 
taking the limit as ri,fi — > —00 and T2,f2 — * 00. 



For parametrization-invariant actions, one can choose to parametrize the trajectory of each particle by Minkowski time, but interaction 
terms involve double integrals over the parameter time of each particle. 
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The variation of the action integral (JSJ) is given by 
W(ri,T 2 ,fi,f 2 ) = mXa -j 4 

/■T2 

+ / drSx 01 

Jti 

+ 



f2 , OA 

CtT-— — 

dx a 



Sx c 



dr (-i 7 iT)^ 



df 



df fe Q 



mi. 



T2 



d <9A 



dr (-x 7 iT)5 



r/r 



8A 



dR a dr dx a 



DA d OA 



dR a df dx a 



where R 



a — a 



R a . 



Requiring that the limit of the variation vanish when Sx a \± oa — and <5a; Q |± 00 = 0, 

lim<5/:= lim SI = 0, 

(Tl,T2)— >(— oo,+oo) 

(Tl,T2)— ► (— 00, +00) 



yields the equation of motion for each particle, 

d mir 



dr (-£ 7 ±t): 



df 



OA d dA 



dR a dr dx a 



(8) 



(9) 



(10) 



mx 



dr (_£ j7) 



r/r 



<9A d <9A 



<9i? a df dx a 



(11) 



The action integral /, with finite values of r^, fj, is invariant under Poincare transformations of the paths that leave 
the path parameters fixed. Invariance of / under the infinitesimal spacetime translation of each path by a constant 
vector a Q , 



Sx a = Sx a = a° 



(12) 



implies conservation of 4-momentum: That is, assuming the equations of motion, (|10|l and i|ll[l . and substituting 
Eq. JHJ) in Eq. ©, we obtain 



51 



(13) 



where 



P a (r,f) 



{—x 1 x 1 ) 2 



df 



dA 



— OO J — OO " T 



(r) 

dA 
~dR* 



mx c 



(— x 7 i 7 ) 2 
drdf. 



, <9A 
dr— — 



(r) 



(14) 



Translation invariance, (5/ = 0, implies that P a is independent of r and r. 

Similarly, invariance of I under an infinitesimal Lorentz transformation of each path, 



Sx a = e^xp and Sx a = e^xp, 



a. ot{3 - 



(15) 



where e a ^ = —t? a is a constant antisymmetric tensor, implies conservation of angular momentum. Again assuming 
the equations of motion, ((THf and |TT|. and substituting Eq. lfT5| in Eq. JBJ), we have 



SI 



= L "/5( T 2' r 2) - Lap(Tl,Ti) 



(16) 
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where 



L*p{t,t) 



m (x a Xp — XfjX a ) 
(— X^X 1 )^ 

ffl (x a X/3 — Xj3X a ) 



, , <9A dA 



dA dA 



OO pT 



OO J —OO J T 



dA 



dA 



Ja dRP X/3 dR a 



Here, Lorenz invariance implies 



e af3 I x p 



dA 
~dR« 



■ X(3 



dA 
dR* 



dA , dA 



dx° 



(r) 
if) 

. dA 

Xa d±p 



= o. 



dA 



drdr 



Finally, 51 = implies that L a p is independent of r and t. 
With a definition 



w := R a R a = {x- x) 2 , 



(17) 



(18) 



(19) 



an interaction A that depends on the positions only through w, A = A(w,x,x), is a restricted form of interactions 
that satisfies the property ©. For such an interaction term, since dA/dR a = 2R a dA/dw, the last line of Eq. [fTT|l is 
written 



OO J — OO J T 



dA 



1 



dA 



drdr. 



(20) 



III. ACTION AT A DISTANCE THEORY FOR POST-MINKOWSKIAN GRAVITY 

As usual in a linearized framework, all tensor indices will be raised and lowered by the flat metric r] a fj of the 
background spacetime. 

A. Fokker actions for point-particles in post-Minkowskian gravity 



Havas and Goldberg [l|| derived equations of motion for point particles in general relativity by expanding the 
metric and demanding that the covariant conservation law for the stress-energy tensor be satisfied to first order in 
the perturbation, with a time-symmetric, half-advanced + half-retarded field for the first-order metric. They found a 
Fokker action, an action integral / for which limW = (in the sense of Eq. ©) gives the the equation of motion to 
the same order. In our notation, their interaction term is 



A(w,x,x) = 2mmd(w) 



(XqX ) — t - ^XqX — |— ^X flX^ ~\~ 2 

(— i 7 i 7 )5(— xsx 6 )^ 



(21) 



Ramond 20] subsequently formulated a general action-at-a-distance theory that is invariant under reparametrization 
and includes tensorial interactions. Following Ramond's argument, we find that the following interaction term is 
reparametrization invariant: 



A(w, x, x) = 2mm S(w) 



(x a x a f 



2 *^ Gt *^ i^iJC 



{— i 7 i 7 ) s (— xsx s ) 2 



(22) 



When parameters r and f are chosen to satisfy x a x a = —1 and x a x a = —1, as in the Havas-Goldberg treatment, 
the equations of motion agree with theirs. Because x a is normalized with respect to r) a /3, n ot r\ a p + h a p, the Havas- 
Goldbcrg r is not an affine parameter. 

We can see as follows how, starting from the post-Minkowski equations of motion, one arrives at an interaction 
term with A given by Eq. (|22J) . The derivation also shows why one obtains a Fokker action, not a true action. 
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To linear order, the metric g a p is a sum of the flat metric, rj a p, an d & perturbation h a p, 

g a p = Vuf3 + h a p. (23) 

where the perturbed metric h a p is a sum of the half-advanced -fhalf-retarded field of each particle. As mentioned 
earlier, in the post-Minkowski treatment of point particles, the effect on the motion of particle m from the field of m 
itself is simply a mass renormalization. That is, the motion of m is described to linear order in h a p by the linearized 
geodesic equation with h a p replaced by the value at the position of m of the half-advanced+half-retarded field h a p 
of m alone: When r is an affine parameter (proper time with respect to the perturbed metric), the geodesic equation 
to linear order in h a p has the form 

ij] a p + h a p)x p + Cc^iPx 1 = 0, where C a p 7 := ~(V,a/i Q7 + V 7 /i /3ct - V a hp 7 ). (24) 



Equivalently, 



^ [{v a p + h^ifi] - ^JvhfriPxi = 0. (25) 



To linear order in h„, the expression (r\ a p + h a p)x^ in Eq. I|24|l can be replaced by x a — r\ a p'x^ ', because x a is already 
order h . The form given in Eq. (|25|l . however, conforms to that of the action integrals below. For t a generic time 
parameter, the geodesic equation has the form given in Eq. I|32|) . This can be directly shown from H25[). but we derive 
it below from the action for a point particle, written to linear order in h a p. 

In the deDonder (harmonic) gauge, V ph a P = 0, h a p is given by 

□ (h a p - ^Vaph) = -16nf a p, (26) 
where the stress-energy tensor T a P of m is defined by 

/oo —a —3 

df8(x-x{?)) XX r . (27) 
-oo ( — X 7 X~t)2 

Here V a is the covariant derivative operator of the flat metric i] a p, and □ = V Q V Q is the corresponding flat 
D'Alembertian. Using the half-retarded + half-advanced Green function, G(x,x) = 6(w) (a solution to DG(x,x) = 
— 4ttS(x — x)), the solution to Eq. (|26|l is written, 

d?6(w)— , 2V f 7 . (28) 

-oo (-X$X S )2 

The trajectory of the second particle, (rh, x a (f)), is similarly a geodesic of a background spacetime with metric 
9aP = Vap + h a p, given by Eq. (|25|l or Ij32(l . with barred and unbarred quantities exchanged. The source for h a p is 
the particle (m, x a (t)): 

/oo +a-p _ l n ap- -7 
dr5(w)— 2 ^ , 7 . (29) 
-oo (-XgX S )2 

(Note: We use a bar to label the perturbed field acting on particle m; h a p is not the trace- reversed form of h a p.) 

To find an interaction term that reproduces the equations of motion, we begin with an action I m for the first particle 
in the field of the second, a linearized geodesic equation on a background spacetime with metric g a p = r\ a p + h a p. 
The action I m is then the action for geodesic motion, 

T2 /-T2 i 

/ An- \-(rr, „-Lh n \T a T^ 2 



— m I dr (—g a px a x 13 ) 2 = — m j dr \_—(rj a p + h a p)x a x 
written to linear order in h a p: 



t i 1 f T ' 2 ± a xP 
I m = in / dT(-i a x a )2 + -m dr h a p- - r . (30) 

2 J T1 (_a; 7X 7)2 



Tl 
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The action I m is invariant under time reparametrization. From its variation, 5I„ 
linearized geodesic equation with arbitrary parameter t: 



0, we can directly compute the 



Sim 



1 



±~<± S 



2 (-x^y. 



■r-0 



(— £ 7 ± 7 )^ 



-Sx a 



dr 



(It 



1 x 7 x s \ x^ 
Vap + Kp + -r) a0 h lS - . —j - . —r\ 

2 (-x^xf)^ J {- x x iy 



1 X^X"! 



6x a . (31) 



Then, requiring SI m = for variations Sx a with fixed end points yields the equation of motion, we have 



d 



I X~^ X^ \ X 1 ^ 
Vaf3 + Kj3 + - VaphiS- . — T ~ ; TTTT 

. 1 y-x^x 1 )^ / (— x-yX 1 )? 



x^x 1 



^ (-.T 7 X 7 )- 



If in the action I m of Eq. (|30|l we substitute Eq. I|28|) , we obtain 

2mm 



T2 i f T2 f°° (x a x a ) 2 

I m = — m / dr (— i Q i Q )2 + 2mm / dr d? 8(w) 



Tl 



i^XqX x^x^ 



(— i 7 iT)3 (-xsx 6 ) 7 - 



(32) 



(33) 



To obtain an action whose variations with respect to both x a and x a yield the equations of motion for each particle, 
one cannot simply add a kinetic term for fh and a second interaction term in which the barred and unbarred quantities 
are interchanged, because the new interaction term would involve x a and alter the equations of motion of the first 
particle. (Roughly speaking, one would be double-counting the gravitational binding energy.) Instead, one observes 
that the interaction term becomes symmetric under interchange of the two particles m, x a (r) and fh,x a (f) if, in the 
infinite integral, we make the replacements — oo — > f\ and oo — > T2- The resulting action integral, symmetric under 
interchange of two particles, has the form of Eq. JSJ, with interaction term l|22|l 



2 ex. X Xj^X 



i f f2 . . i f T2 f^ 2 (x a x a ) 2 

dr (— i Q i Q )2 — fh I df (— x a x a ) 5 + 2mm / dr I d? 5{w) — 

Jfi Jn J fx (— i 7 iT) 2 (— xsx 8 ) 



(34) 



The price for this symmetry is that a variation of I yields equations of motion in which the integrals giving h a and 
h a /3 extend only from fi to T2 (or n to t%). It is for this reason that action-at-a-distance theories require a Fokker 
action, whose equations of motion are given by \m\5I = 0. 

As mentioned, the above theory is invariant under reparametrization of r and f. The obvious choice of a proper-time 
parametrization, with r\ a ^x a x^ — —1 and fj a ^x a x^ = —1, simplifies computations of physical quantities. 



B. A Fokker action for the affinely parametrized equations 



One can, of course, specialize the parametrization-invariant action to affinely parametrized trajectories, but afhne 
parametrization also allows a generalization to interacting particles of the quadratic action for geodesic motion, 
IT 2 .. 

— m / drx a x a . To construct an action integral / with quadratic kinetic term, for which lim<5/ = reproduces the 
2 J Tl 

affinely parametrized equation of motion, Eq. (|25|l . and its barred «-> unbarred form, we modify the kinetic terms in 
the action integral of Eq. (JSJ , writing 



I{ T i, T 2, ti, T2) = -m / dr x a x a + —fh / df x a x a + / dr d?A(w,x,x); 

2 J Tl J T\ Jt\ J fx 

and we take as the interaction term 

A(w,x,x) = 2mfhS(w) 

Affine parametrization, the requirement 

(Vaf3 + h a fi)x a xP = constant, 



{x a X ) ^XqlX XfiX r 



(Vaf3 + h a p)x a x = constant, 



(35) 



(36) 



(37) 



is enforced by the equations of motion. 



8 



Formulas for the equation of motion Eq. (|10(l and Eq. <|11|) . the 4- momentum Eq. (|14fl and the angular momentum 
tensor Eq. I|17l) are changed as a result of this modification of the kinetic terms in Eq. (|35[) . The variation of our 
kinetic terms 



IT 2 .. 

— m8 j dr x a x a = mx a 5x° lT 



Tl 
T2 



dr mx a 5x a , 



—m8 I dr x a x a = mx a 8x c 



Tl 
T2 



dr mx a Sx a , 



(38) 
(39) 



results in the following form for the equation of motion, 

OA d dA 



rax 



df 



dR a dr dx° 



(40) 



dr 



dA d dA 



dR a df dx a 



~(h al3 x p ) + 'ahp-fX^ 1 x 1 ', 



in agreement with the linearized geodesic equation lj25|) . 
The 4- momentum is given by 



P a (r,f) 



and the angular momentum by 



m:r 



df 



dA 



(r) 



, dA 
dr—r— 

dx a 



OO /» T 



dA 

dR* 



drdr, 



m (x a ±f3 — xpx a ) + 

m {x a X{3 - XpXa) + 



dr | .c, 
dr I x, 



dA dA 

a diP~ X0 dx^ 
dA _ dA 

d¥~ Xp d¥ 



(r) 



OO p T 



— OC J — OO J T 



M( xp x a -x a x )-\\± fi j^-.r.,— , 



dA dA \ 

)- Xr.~ 



drdr, 



(41) 



(42) 



(43) 



These expressions for momentum and angular momentum are formally identical to Eqs. I|14|) and 11711 if one sets 
{—r] a px a x 13 ) to 1. This is misleading: In the present section, x a is a unit vector not of the Minkowski metric but 
of the perturbed metric r} a p + h a p, and the altered form of A compensates for the altered normalization, leading to 
identical Newtonian limits of the two expressions. 



C. Circular orbits 



We consider now a system of two point particles in circular orbit, for which the gravitational field seen by each 
particle is the linearized half-advanced + half-retarded field of the other. This is a gravitational analog of a solution 
obtained by Schild |l If for two point charges. In the gravitational context, however, linearized gravity is a toy model 
for bound systems: Discarded nonlinear terms are of the same post-Newtonian order as linear terms that are kept, 
and hence of the same magnitude as terms of the next post-Minkowskian order. 

We introduce a basis {t a , w a , <fi a , z a } of unit vectors of the flat metric r) a p. A particle in circular orbit in the z = 
plane with constant orbital radius a has cylindrical coordinates {t, zu = a, 4>, z = 0} and position vector 

x a =tt a +am a . (44) 

Its spacetime trajectory is tangent to the helical Killing vector 



(45) 
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where (f> a — m<f) a is a rotational Killing vector of the flat metric r) a p and f2 is the particle's constant angular velocity. 
With 7 := dt/dr and v :— ail, the particle's 4-velocity and acceleration are given by 



and 



x a = 7 fc Q = j(t a +V(j) a ) 



The second particle m has a circular orbit of radius a about the same origin with coordinates 
{t, w = a, <p, z = 0} and position vector 

x = tt + avo . 

The particle trajectory is again tangent to the helical Killing vector 

k a = t a + nj) a , 



(46) 
(47) 

(48) 
(49) 



with 4> a the value of the vector field <j) a at the position of the second particle. With 7 = dt/dr and v := afl, the 
acceleration and 4-velocity of the second particle are given by 



and 



x a = jk a = j(t a +v(f) a ) 



x a = -fjjfiw Q 



(50) 



(51) 



For the parametrization-invariant formulation, it is convenient to parametrize the trajectories by proper time with 
respect to the flat metric rj a p] their tangent vectors are then unit vectors of 7? Q /3, 



r] al3 x a x 13 = -1, T] al 3X a x fJ = -1. 



(52) 



The quantities 7, v, 7 and v obey the familiar relations 

7=(l-u 2 ri 7=(l-u 2 )-3. (53) 

For the affinely parametrized formulation, a correct Newtonian limit is reproduced when a normalization with the 
perturbed metric (Sec. IIIIBf) is used: 

(Vo C p + h af3 )x a x = -1, (fj a p +h a p) x a x = -1, (54) 

which yield relations 

7 = (1 - v 2 - h^k 13 )-^, 7 = (1 - v 2 - h aP k a ¥)-^. (55) 
With the first particle, m, at <j> — for t — 0, its trajectory has coordinates 

<j) = Qt , t = 7 r; (56) 

and the trajectory of m then has coordinates 

4> = 7T + Cli , i = 7f-7r/0. (57) 

The positions of both particles are specified by a single parameter, and it is natural to choose a descriptor r\ of their 
relative motion, defined by 



where we pick r = r = when <f> = <p = 0. 
A vector R a :— x a — x a becomes 



R a 



1 

a 
1 



rj := (p — <j) = 17(7T — 7r), 

(7r — rf) t a + (v — v cos rj) w a — v sin rj 4> a 
(ir — i])t a — (v — v cos rj) m a — v sin rj (j> a 



(58) 

(59) 
(60) 
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The half-retarded + half-advanced Green function S(w) = 6 [(x — x) 2 ~\ has support on the two events that correspond 
to the roots of 

10(77) : = i x ~ x) 2 = ^2 ~ + v 2 +v 2 - 2iwcosr/] = 0. (61) 

These roots are given by 

T] = 7T ± ip, (62) 

with if the positive root of 

ip 2 = v 2 + v 2 + 2vv cos ip. (63) 

Here ip + 7r is the angle between m and the retarded position of m, and Eq. i|63|) has a simple geometrical meaning, 
illustrated in Fig. 

95 + 7T = rot - (/>, 

in which R is the distance from m to the the retarded position of m. That is, using the geometric relation 




FIG. 1: The figure shows mass m at time t and mass m at the corresponding retarded time. The angle re t is the <f> coordinate 
of the intersection of the past light cone of m with the trajectory of m, when m is at <j>. 



R 2 = a 2 + a 2 - 2aacos(4> Ict - 4>) = 7^7 (w 2 + v 2 + 2vv cos tp) , (64) 
noting that R = i lct — t, and using Eqs. (|56I57() . we have 

<P = 0rot —7T—(p = fi(?ret ~ *) = ^-R- (65) 

Then (p 2 = n 2 R 2 and Eq. I|63|) is immediate from l|64|). Analogously, the root r\ = n — <p corresponds to the advanced 
position of m defined by — tp + n = a d v — <j>- 



D. Conserved quantities and angular velocity of circular orbits 

We now rewrite the integrals appearing in Eqs. QlOfl. I jllj l. (|14|l . and l|17(l for a parametrization-invariant action, 
in terms of the parameter 77. To evaluate the integrals, we use formulas obtained in Appendix IA II The resulting 
formulas for affinely parametrized formulations Eqs. 1)41)0. (|41|1 . I|42l) . and (|43|l are quite similar to those for the 
parametrization-invariant formulations. However, the fact that the action appropriate to an initial choice of afhne 
parameters is not a special case of the parametrization-invariant action leads to definitions of a conserved energy that 
are not identical in the two formulations. 
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1. Parametrization-invariant model 



For circular orbits, only the radial component of the equations of motion is nontrivial, and it has the form, 

7 OX/3 



-iniiir/-. j dr) — j 5'(w)(v — v cosr/) $(77, v, v) 



-S(w) [(1 — vv cos rf)v cos 77 — — v(l — v 2 ) 



1 v 



2 1 -v 



m^vfl = 



A I dA 2 , 



7 dA 

7 3£/3 



(66) 

(67) 
(68) 



= 477777177 / dr] < — 6'(w)(v — v cos 77) $(77, u, v) 



-S(w) [(1 — vv cos 77)77 cos 77 — — v(l — v 2 ) 



1 V 



2 1 - v- 



;$(T],V,V)]\, (69) 



where Eqs. (|A10|I and (|A1 1|) in Appendix IA II are used. We have introduced here a function <b(rj,v,v), defined by 

$(r?, v, v) := \(x a x a ) 2 - \x a x a Zp&] = (1 - vv cost?) 2 - 1(1 - t; 2 )(1 - tj 2 ). (70) 



Since the center-of-mass frame is chosen, the only nonzero components of the 4-momentum and the angular mo- 
mentum are E := —P a (T,?)t a = —Po(t,t) and L := Li2(t, f). Taking (r, f) = (0,0), we have 



£ = P°(0,0) 



1 



<9A 



• -7T / -77- 

7" 0^0 



(0) 



TflX 



1 

7^ 



r/r/ 



o»A 

9i 



(0) 



1 



77f2 2 



d77 77^2(2;° -5°) 



in", FT~^" / drjS(w) [(1 — vvcosrj) — — (1 — t; 2 ) — — j$(r7, u, w)] 



7717 



n 

477777777 
477777777 



L = L 12 (0,0) 



m(xiX2 — x 2 ii) + 



m(xiX2 — X2X1) + 



dr]S(w)[(l - vv cos 77) - 1(1 - v 2 ) - _ 2 $(r),v,v)] 



drj 5' (w) 77(77 — 77) $(77, t>, v), 

OA dA 



(71) 



1 

7T2 
1 

7~f2 



, . dA dA 
dr] X-L—zt - X2-T- 
x dx z dx L 



(0) 
(0) 



77SI 2 



^77 77 



dA 

dw 



rwyv Ammjjvv 
rfir^v 2 Amfhjjvv 



If dA dA\ 

{x 2 X 1 -X 1 X 2 )--[x 2Wl -X 1 —j 

1 v 

dr\ 5(w) \(1 — vv cos 77) cos 77 — — — (1 — v 2 ) 

2v 



(72) 



lv 1 



2 v 1 — v- 



-$(r7, v,v)] 



n 2 

4,mfh"f'yvv 

02 



d-q S(w) \(1 — vv cos 77) cos 77 — — — (1 — v 2 ) — — — $(??, v, v)] 

2 v 2 v 1 — ir 



d>] 



O 2 



5'(w) 77 sin?7 $(77, v, v) + 8(w)(l — vv cos 77)77 sin 77 



(73) 



where we have used a formula, Eq. I|A5(1 . derived in Appendix IA II 

Formulas for integrations involving a 5-function are summarized in Appendix IA 21 First we define a function 
$(95, v, v), related to $ by an integration: 



1 ,' - t 1 : < < v v I,; „ ; 77 

&{<p,v,v) = I dr]5{w) 9[Tj,v,v) - 



(1 + zwcos^) 2 - i(l - w 2 )(l - v 2 ) 



if + vv sin ip 



(74) 
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Using this relation, we can evaluate the integrals appearing in the equations of motion l|67|l and l|69|l . obtaining an 
algebraic relation for the radius of each orbit (or, equivalently, for the velocities v — ail and v — att) in terms of the 
angular velocity Q: 

— my 2 vfl = 4mm7 2 7f2 2 ' | (1 + ww cos <p)v(<p> cos w — v 2 sin w) H — w(l — v 2 )((p + ww sines) 

(ip + ww simp) z i 1 2 

1 v 1 

\v simphp + vv simp) + (1 + vv cos ip)(v + v cos <y9) -(tp + vv sin<^) 2 l $(</5, w, w) > , (75) 

2 L ' \ — v A J J 

— 7f77 2 wf2 = 4mr7i77 2 r2 2 - — ; I (1 + vv cos (p)v(ip cos ip — v 2 sin ip) -\ — w(l — v 2 )(<p + vv simp) 

[}P + vv simpY y ' 2 

1 v o 1 

I"w sin w(w + sinw) + (1 + «5 cos + v cos <p) =z(<P + vv siniyj) 2 ] §Up, w, w) > . (76) 

2 L ' \ — v 1 J 

The expression for the angular momentum l|73[) is integrated by substituting the equations of motion Eq. (|67|l and 
(EU into (J73J, 



477777777WW 
2 



Explicitly, 



dr/ <5(u>) (1 — vv cos 77)77 sin 77 

dr\ 5'(w) $(77, w, w) (w 2 + w 2 — 2vv cos 77 — vvr] sin 77) (77) 

o 

d?7 [2wS'(w) + 5(w)] $(77, v, v) (78) 
dr/ <J(w) $(77, t>, v) — 2mfh / jj $>(<p, v 1 v). (79) 

L = 27,77777 (1 + ^ COS ^^ (1 ^ 2)(1 '' 2) , (80) 
(/? + ww sin </? 



_4mm77 y/ 00 JmV ^„^^ „, _ W/i>2 , ... 
2mfh r fj 

277777777 

n 2 



with given by Eq. 1)63(1. 

To compute the energy £7, we first compute 75 — f2L: 

E-VtL = 7777(1 -w 2 ) + 7777(1 -w 2 ) 

477777777 Z" 00 - „ 

— — / drj8(w) $(77, w, w) + (1 — wwcos?7)ww77sin?7 

477777777 /*°° ~ 

— — — - — / drj 6'(w) 77 [(7r — 77) + wwsin?7] $(77, w, w) (81) 

^ ' J — 00 

= 7777(1 - w 2 ) + 7777(1 - v 2 ) - 2m ™ 7 ^ f dr]S(w) $(77, w,w) (82) 
= 7777(1 — w 2 ) + 7777(1 — w 2 ) — 2mfh'y"ffl w, w). (83) 
Using the definitions (|53|l of 7 and 7, the expression for E takes the simple form, 

„ 777 777 , o, 1 , n. 1 . 

E =— + — =m(l -w 2 ) 2 +777(1 -w 2 )2. (84) 

7 7 

As noted in the introduction, this expression for E is identical to the expression for E in the electromagnetic two-body 
case derived by Schild [llj, and to that for a scalar field, obtained in Appendix[U] To clarify the structure of relevant 
equations used in calculations above, their formal expressions are presented in Appendix ID 21 

The corresponding expressions for angular momentum are each proportional to the potential dotted on each free 
index with the helical symmetry vector k a . For a scalar interaction, masses m and 777 have scalar charges q and q. A 
scalar field ip at m due to q satisfies 

Difj = -4nq j dfS 4 [x - x(f)], (85) 
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with the corresponding definition of tp. Similarly, for an electromagnetic interaction with charges e and e, the vector 
potential A a at m due to e satisfies (in the Lorentz gauge) 

UA a = -4ne f d? u a (?) 5 4 [x - x(?)]. (86) 

With ip, A a and h a p evaluated at x, if), A a and h a p evaluated at x, the angular momentum L has the following forms 
for scalar, electromagnetic, and post-Minkowskian gravitational interactions: 

Scalar charges q and q 
Electromagnetic charges e and -e 
Post-Minkowski masses m and m 



L 
L 
L 



1 i 
--if, 

7 



7 



q 2 1 

77 tp + vv cos tp 

1 + VV COS C2 



(87) 



— —A h a —~—A h a — pp 

Si Si </? + VV COS </3 



2Q 



2sl 



a/3* 



(89) 



2. Equations of motion, energy and angular momentum for the affinely parametrized model 

The analogous computations for the affinely parametrized formulation are outlined here. Performing the integrals 
in the radial components of the equations of motions (|40(l and (|41|) . we obtain 

— mj 2 vD, = 4mfhj 2 jCt 2 < (1 + vv cos <p)v(<p cos tp — v 2 sin p) H — v(l — v 2 )(tp + vv sin tp) 

(tp + to ship) 2 ^ 2 

— — [w siny(y + sin tp) + (1 + wu cos y)(i> + w cos y)] $(co, u,£>)l , (90) 

— fhn 2 v£l = 4mfh'y^ 2 VL 2 — { (1 + vv cos p)v(pcos tp — v 2 sin tp) + —v(l — v 2 )(tp + vv since) 

(tp + vv sm tp) z y ' 2 

— — [v sin tp(tp + vv sin tp) + (1 + vv cos tp)(v + v cos co)] $(co, j . (91) 

The angular momentum turns out to have the same form as in the parametrization-invariant formulation, 

_ (1 + vvcosp) 2 - i(l - v 2 )(l ~ v 2 ) 

L = 2mm"f"f^(tp,v,v) — 277777777 ■ (92) 

tp + vv sm tp 

We find, however, that the form of E — f2L differs from the corresponding Eq. (|83|l : 

E — QL = 7717(1 — v 2 ) + 7777(1 — v 2 ) — 67777777751 <&{tp, v, v). (93) 

The normalizations l|54|l have the explicit form 

(r)af) + hafi) ± a x p = - 7 2 (1 - v 2 ) + 4fhj 2 jn ${tp, v,v) = -l, (94) 
(fj a p + hap) x a xP = -7 2 (1 - v 2 ) + AmjfQ <f>(tp, v, v) = -1. (95) 

Substituting these in Eq. (|93|l . the expression for the energy becomes 

E = h — + 47777777711 $(tp, v, v) = 7777(1 — v 2 ) + 7777(1 — v 2 ) — 47777777717 <&(tp, v, v). (96) 

7 7' 

As implied by Eqs. i|92[l and i|93|) . the form of E differs from the energy l|84[) of the parametrization-invariant action. 
Both, of course, agree in the Newtonian limit. 
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IV. FIRST LAW OF THERMODYNAMICS FOR BINARIES DESCRIBED BY FOKKER ACTIONS 



The first law of thermodynamics governs nearby equilibria of conservative systems. For binary systems with circular 
orbits, an equilibrium solution is a solution that is stationary in a rotating frame, a solution with a helical symmetry 
vector k a . When internal degrees of freedom (e.g., baryon number, entropy, vorticity) are fixed, the first law relates 
the change of energy to the change of angular momentum in the manner |5j 



se = nsi. 



(97) 



In the present context, however, the presence of a radiation field whose energy is infinite makes the relation suspect; 
and the lack of a true action means that the simple Hamiltonian proof for point particles with a Newtonian potential 
does not hold. But the relation i|97[l is true, and its proof is an extension of the Hamiltonian proof that uses the 
parametrization invariance of the relativistic Fokker action. 

Proposition. Consider a parametrization- invariant Fokker action of the form © . Suppose there is a family of solutions 
x a (s, r), x a (s, f) for which the particles move in circular orbits with angular velocity f2(s), with spacetime trajectories 

along the helical vector field k a (s) =t a + tt(s)(j) a . Then 6E = QSL, where SQ ■ ' K * 



ds 



s=0 



Parametrization invariance and constraints. 

As a prerequisite to the proof, we begin with a brief discussion of parametrization invariance of an Fokker action 
of the form J5J and the constraints that follow from it. These are analogs of the Hamiltonian constraint (vanishing 
of the superhamiltonian) associated with a true parametrization-invariant action. A parametrization r — > /(t) maps 
the path x a (r) to the path X a — x a o /. Similarly, f — > /(f) maps x a to X a — x a o /. Invariance of the action, 



I=- dr m(-x a x a ) 1/2 - df rh(-x a x a ) 1/2 + 

under reparametrization follows from the scaling 

A(X — X, X a , X a 

For /(t) = fcr, /(f) = f , we have 



dr 



dr A(x — x, x, x), 



/(to)/(t ) A{x - x,x a ,i a )\ /(ro) , /(fo) 



(98) 



A(X - X,X a ,X a ) 



To/fc, T 



A(x ~ x, kx a , x a )\ ro ^ o = k A(x - x, x a , x a )\ ro ^ fo . 



Then the relation, 



— [kA(x — x, x, x)] 
dk 



k=l 



— [A(x — x, kx, x)] 
dk 



k=l 



and its barred <-*■ unbarred counterpart imply 

A = x a 

Proof of Proposition. 

The equations of motion involve the potentials 



dA 



A = 



, dA 
dx*' 



(99) 



Un 



drA, 



drA. 



That is, with 1-particle Lagrangians C m and Cm defined by 

£"m v?'i( x x a *} ^ -\- JA mi Cffi m( x x a ) ^ -\-lAfii 
Eqs. I |l Up and (|ll|l are the equations of motion for the actions 



drC r , 



drCf, 



(100) 



(101) 



(102) 
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' * 0, _ ™ -_2= . (103) 



namely, 

dr V dx a J dx a ~ df \ dx a J 3x a 

The one-particle 4-momentum p a associated with the action I m is 

_ dC m _ mx a dlA m (-\r\A\ 

Pa ~~di^~ (-iT± 7 ) l /2 + a±5"- liU4j 

Because £ m and i" have the same reparamctrization scaling, namely 

C m ^fC m , x a ^fx a , (105) 

the momentum p a is independent of the choice of parameter. 

Let t be a choice of Minkowski time, t a the corresponding Killing vector (<9 t ), and cf> a a rotational Killing vector 
orthogonal to t a . The 1-particle energies associated with I m and I m are 

<9£ m mi 9W m 

^ = - tPa= -^r = F^)^-^r ; (107) 

and the 1-particle angular momenta are 

L m = <j> a p a = L fi = 4Tp a = (108) 

We first use the scaling relation and the equations of motion to show that these 1-particle momenta and angular 
momenta satisfy SE m — QSL mi 5E m = £l8L m . The scaling of C m and C m (parametrization invariance of I m ) implies 

(In other words, the 1-particle superhamiltonians vanish: Tt m := x a p a — C m — 0, 7i m '■— % pa — C m — 0.) Then 

x a BC C BC C 

t 3x a t 3x a t v ' 

(Note that C m /t is the form of the Lagrangian appropriate to Minkowski time: J dtC m /t = J drC m .) Consider now a 
family of solutions to the equation of motion for m with circular orbits, each solution stationary in a comoving frame. 
Neighboring orbits satisfy 

BC BC BC BC 1 BC 

where the equilibrium condition — was used to infer Sx a ^f^- = 0. 

era 7 ® x 
The total energy E of Eq. (|14|l is not, however, the sum of the 1-particle energies: E m + E m does not include the 

field energy and is not, in general, conserved. Instead, the total energy has the form 



E = -P a t a = E m + E m - I dr l"Y n U n , - I ,/rrV,,//,,,. (112) 
The total angular momentum is similarly 

L = P ( A a =L m + L m + I dr <j) a V a Um + I df (j) a V a U m , (113) 



where we have used the notation 

f _ BU m _ 

dx a 34) ' 3x a di> ~ ~ T Bx a ' r 3x a 
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(See also Eq. lfD23)> of Appendix El) 

To recover the first law, one uses the fact that E and L are independent of t and f , setting r and f to — oo to 
eliminate the final term: 

SE - flSL = lim (8E - Q8L) 



T.T — >* — OO 



= lim 



T.T — > — 00 



SE m - £l5L m + SE 7 - n - flSL m + / dTk a \7 a 5U m + / dfk a V a SU 7 - n 

J — oc J —oo 

= 0. □ (115) 
Note that the Killing symmetry, k a V a Um — 0, does not in itself imply k a V a 8U m = 0, because Sk a — (f) a 8Vl ^ 0. 

First law for afRnely parametrized action 

The first law also holds for the affinely parametrized action i|35f) . The proof is similar to that for the parametrization- 
invariant action, with one-particle Lagrangians now defined by 

in m 

with the affinely parametrized A replacing the parametrization-invariant A in the definition of Li m and him ■ Note that 
affine parametrization is equivalent to the conditions 

2C m = -m, 2Cm = -to, (117) 

while the fact that C m and C m are quadratic in the velocities x and x implies 

2C m = Pa x a , 2C m =p a x a . (118) 

The first of these relations, (|117|) . together with the equilibrium condition dC m /dvj = 0, implies 

5£ m = ^Sx a = Pa Sx a . (119) 

From the second relation, l|118|) . we have = —8m = S(p a x a ) = Sp a x a + p a Sx a , whence, by (|119fl . 

x a Sp a = 0. (120) 

Again writing E m = —pt, x a = m'w , we have 

SE m = v a 5 Pa = toSLm. (121) 

Finally, the relation 

5E = Q5L (122) 
follows as before from Eqs. I|112|) - (|115|l . which hold as written. 

V. ACTION ACCURATE TO FIRST POST-NEWTONIAN ORDER 

In a first post-Newtonian approximation, the equations of motion have corrections of order v 2 /c 2 to their Newtonian 
terms. As wc noted above, the first post-Minkowski approximation, by including only terms linear in m and m, fails 
to be accurate to first post-Newtonian order. The omitted terms in the equation of motion are quadratic in the masses 
m and m, and they arise from a single term in the post-Newtonian Lagrangian, namely 

mfhim + m) 

LpN = (123) 

We can make the Fokker action accurate to first post-Newtonian order by adding any term that agrees to this order 
with H128JI . and we have tried two alternatives: A simplest parametrization-invariant choice of Apn that reproduces 
CM is 

A PN = -£ (i-t) K — 2 Hx^xsx 6 ) 1 / 2 , (124) 
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with value for a circular orbit given by 



fl 3 mm (to + to) 
2 v 2 + v 2 — 2vv cos rj 



~ , , u u It bl lb \ lib II b I 

A PN = -6 (r, - tt) — a| „ - v (125) 



A choice that is special-relativistically covariant (and parametrization invariant) is 

a xi ~,mm(m + m) (x^x^xsx 5 ) 3 ^ 

AsPN = 2 {Ra±aRfi &y,2 - ( 126 ) 

The corresponding corrections epN and esPN to the conserved energy and the corresponding corrections ^pn and 
^SPN to the angular momentum are given by 

e PN = ^OIp N , (127) 



mm(m + m)n 

*pn = —, — _s 2 , (128) 

77(w + w) z 



e S PN = ^^SPN, (129) 
mm(m + m)Q 1 

(130) 



t-SPN = 



(^)3/2 (Lp + vvsmip) 2 

The derivation of Eqs. l|127|) - (|13Ufl and corrections to the equations of motion are given in Appendix iBl 



VI. NUMERICAL SOLUTION OF ORBITAL EQUATIONS 



We have numerically solved the orbital equations, Eqs. 17511 and l|76(l for the parametrization-invariant model, and 
Eqs. I|90|l and (|91|l for the affmely parametrized model, finding a(f2) for equal-mass particles. Figs. |21 and [21 show the 
relations for the post-Minkowski action. 



— PM inv. 


i-C/271 


- - PM affine 


r=C/27i 


— - PM+SPN 


r=C/27i 


PM inv. 


r=a 




1U 0.0 0.1 0.2 0.3 0.4 

V 

FIG. 2: The velocity v of a circular orbit versus the radius r, written in the dimensionless form r/M, for particles of equal mass 
M := ra + m, where the r is either the proper circumferential radius C/2tt or the radial parameter a. The motion is governed 
by the uncorrected parametrization-invariant post-Minkowski action (PM inv.), the affine parametrized action (PM afHne) or 
the parametrization-invariant post-Minkowski action with a covariant post-Newtonian correction term 112611 (PM+SPN). For 
the parametrization-invariant case both C/2tt and a are shown. 

As v increases, relativistic beaming decreases the strength of the gravitational field due to m at the position of to, 
and vice- versa. In the case of scalar and electromagnetically bound charges, the smaller held leads to a sharply smaller 
radius for a given velocity v; and the same effect implies a larger value of fi at fixed v. Gravity, however, has the 
competing effect that, at small radius the field is stronger than the Newtonian field. In the exact theory, the result is 
that, for circular orbits about a fixed mass, the relativistic relation between f2, r (and a velocity defined by v = fir) is 
identical to the Newtonian relation, when r is taken to be the circumferential radius. In the post-Minkowski models, 
as Figs.0and|31 the outcome depends on which post-Minkowski action one chooses. For the parametrization-invariant 
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0.0 0.1 0.2 0.3 0.4 

V 



FIG. 3: The corresponding relation between Q and v for the orbits of Fig. [21 with fl written in the dimensionless form Q.M. The 
dashed curve (PM+SPN) shows the small correction to the solid curve (PM inv.) that arises from adding the post-Newtonian 
correction term l|126|l to the post-Minkowski action. 



1.0 




QM 

FIG. 4: The corresponding relation between Q and \h a px a x l3 \ for the orbits of Fig.|5] The dashed curve (PM+SPN) shows the 
first order post-Newtonian correction that is nonlinear in m 1126H weakens the gravitational field. 

action, relativistic beaming dominates, giving a smaller radius for the same value of v and a correspondingly larger 
value of at fixed v. The affine action, in contrast, gives a larger value of r at fixed v, an effect so pronounced for 
v > 0.3 that r reaches a minimum value and then increases with increasing v. 

We were also surprised by the fact that the nonlinear term in m and m that is absent from our first post-Minkowski 
action has sign opposite to the Newtonian potential term. Fig. |3|and Fig. 0] show only a small correction due to the 
post-Newtonian term (denoted PM+SPN in the graph), and the correction weakens the gravitational field. 

Figures Eland El display the values of the angular momentum and energy as functions of S2, for the parametrization- 
invariant post-Minkowski (PM inv.) action, the affine parametrized post-Minkowski action (PM affine), the invariant 
post-Minkowski action with the special relativistically covariant first post-Newtonian correction l|12tj[l (PM+SPN), and 
the invariant post-Minkowski action with first post-Newtonian correction (|124fl (PM+PN). The values of Newtonian to 
third post-Newtonian approximations (OPN - 3PN) are also plotted for references. Again relativistic beaming appears 
to dominate the post-Newtonian correction in the PM+PN and PM+SPN action, leading to a graph in which the 
energy and angular momenta have no minima. Because the relation dE = QdL can be used to show that the innermost 
stable circular orbit (1SCO) is a minimum of E and L, the ISCO present in the post-Newtonian approximation does 
not appear in PM inv., PM+SPN or PM+PN. Interestingly, PM affine model has a simultaneous minima in E(il) 
and L(fl) for a sequence of circular solutions. Because the relation dE — ildL also holds for PM affine model as 
shown in Sect. lIVI the circular solutions with SIM larger than the value at the minima, QM = 0.0522, are likely to be 
dynamically unstable. The values at the turning point of the other quantities are v = 0.184, E /M — 0.988744, and 
L/M 2 — 1.0279123. Each number above is given to computational accuracy. 

Verification of our numerical solutions can be seen from Fig. \7\ in which 1 — dE/QdL is plotted against a measure 
A of the resolution of dL. As shown in Sect. IIVI the relation dE = UdL is exact for Fokker actions of the form of 
PM inv., PM affine PM+SPN and PM+PN, and the numerical determination finds no discrepancy. (The result of 
PM+PN is omitted in Fig. [7] to avoid redundancy.) 
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'6.0 0.1 0.2 ~ 0.3 
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FIG. 5: Angular momentum, in dimensionless form L/M 2 , is plotted against angular velocity for 8 cases. In the key labeling 
the curves, 0PN-3PN refer to models of equal-mass point-particles in circular orbit in the Newtonian, first post-Newtonian, 
second post-Newtonian, and third post-Newtonian approximations from the bottom to top. 




FIG. 6: Energy, in dimensionless form E/M, is plotted against angular velocity, for the models of Fig. 




FIG. 7: Accuracy of the relation dE — QdL. Thin solid lines have an inclination proportional to 0(A 4 ), where A is proportional 
to dL. dE/QdL is evaluated at v — 0.15 using the 4th-order accurate finite difference formula (Lagrange formula). 



VII. DISCUSSION 



In the previous section, we discussed numerical results for the special case of gravitationally interacting particles in 
circular orbits. Here we consider a few surprises that arise in our more general study of Fokker actions. 
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Conserved quantities for particles with scalar, electromagnetic or gravitational interactions are ordinarily written 
as an integral over the field and its first time derivative on a spacelike hypersurface, together with sum of terms 
involving the position and velocity of each particle at its position on the hypersurface. In a Fokker formalism, the 
integral over the field is replaced by a sum of integrals over the trajectory of each particle. A striking feature of 
the resulting conserved 4-momcntum and angular momentum is that they break up into a sum of quantities that 
are separately conserved for each particle, as in Eqs. I)112|l and (|113|) (or IjDISJl and l|D20J) ). In these equations, the 
separate conservation of each contribution follows from the fact that the contribution associated with each particle 
depends only on the proper time of that particle. How is this possible, when the field energy measures the interaction 
between the fields produced by each particle; how is it consistent with the fact that the total 4-momentum and angular 
momentum are, in general, an exhaustive set of integrals of the motion? 

The answer is that the individual integrals that appear in the sum cannot be written as integrals over a hypersurface 
of a density locally constructed from the field and its first time derivative. They therefore do not represent new integrals 
of the motion. In fact, their existence is really a trivial consequence of the equation of motion satisfied by each particle: 



3£ m 
dr V dx a 



dx a 



= 



(131) 



where as before 



-( 



dfA. 



(132) 



Integrating the equation of motion from t\ to T2, we have 



T2 , dU m 
dr— — = 0, 

dx a 



Pa 



dx a 



(133) 



implying P ma is independent of r, where 



dr 



~dx~° 



(134) 



More generally, dotting a vector into the equation of motion (|131fl yields the identity 



d_ 

rf7 



mx a ( a dU m 
(-ip±P)i dx<* 



c 01 

9U m + £ a 9Um 



dx a 



dx° 



mx a C, c 



{-x^x-r) 1 /' 



= 0. 



(135) 



Then, because a Killing vector £ Q of Minkowski space satisfies x a ( a — x a x^V a C,p = 0, the quantity Q m defined by 



Qm — 



mx a C a ra dU, 

i ~r C, 



(-xpxPy 



dx a 



(r)-l dr[C 9 ^+C m 



dx a 



dx a 



(136) 



is independent of r. In particular, for C, a — —t a and £ Q = <p a , Q m is the one-particle contribution to the energy and 
angular momentum, respectively. The generalization of these relations to an rt-particle action is immediate. 

The conserved quantities that arise from a Fokker action involve integrals along the particle paths whose integrands 
are not perfect time derivatives: That is, the integrals are path dependent. 2 Although the total 4-momentum and 
angular momentum are also expressed as path-dependent integrals in the Fokker formalism, they, presumably, differ 
from their 1-particle constituents by the fact that they can be expressed as integrals over a hypersurface of a local 
function of the field and its first time derivative (together with the kinematic momentum of each particle) . The parts 



2 One can similarly construct a path dependent conserved "momentum" for a single Newtonian particle moving in any external field by 
writing 

dp/dt = -VU(r) 

dVSJU = constant. 

Because the integral depends on the path r(t) from to to t, one does not have a first integral of the equation of motion, in the sense of 
a conserved quantity that depends only on the position and velocity of the particle. 
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of these conserved quantities that directly involve the fields - the sum of terms involving lA m for each m - is not, 
however, the full field momentum or angular momentum, 

J TfCdSp. (137) 

Instead, only contributions to the field momentum (or angular momentum) that arise from the product of fields from 
different particles are present in the Fokker momentum (angular momentum). Terms quadratic in the field of a single 
particle must be absorbed in the renormalized mass. 

This analysis leads to a conjecture for why the conserved Fokker momentum and angular momentum are finite for 
bound systems. When terms quadratic in the fields of individual particles are subtracted, only a finite remainder 
survives. 

Finally, we comment on ambiguities associated with the choice of gauge and with the choice of post-Newtonian 
correction terms. Solutions for two-particle circular orbits presented in this paper are obtained in the deDonder 
(harmonic) gauge. A sequence of equal mass solutions terminates near v ~ 0.40 for the parametrization invariant 
action, and near v ~ 0.32 for the affine action, because in each case the equations of motion have no solution for v 
larger than a critical value v. When post-Newtonian corrections are added to the parametrization-invariant action, 
we find solutions only up to v ~ 0.19. When one proceeds, for example, to higher order calculation aiming to identify 
an ISCO, which may be expected to appear near these terminal values, a choice of gauge may be reconsidered with 
care. 

A post-Minkowski action corrected by the terms Aspn or Apn, does not exactly reproduce the first post-Newtonian 
equations of motion, because the post-Minkowski action contains terms of all post-Newtonian orders in the velocities. 
A more systematic approximation to the higher order post-Minkowski Fokker action may be obtained by deriving a 
form of interaction term that agrees with a formal expansion of the stress-energy tensor (see, for example, plf'). 
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APPENDIX A: FORMULAS FOR INTEGRATION 



We present here a formalism for evaluating the integrals that arise in the equations of motion, and in the expressions 
for momentum, energy, and angular momentum of point particles in circular orbits. This is closely patterned on 
formulas used by Schild in [22j in the case of two electrically charged particles. 

1. Change of integration variables 

We first derive a relation for rewriting integrals with respect to proper time r and f as integrals with respect to 
the parameter -q. From Eq. iJSS}, 

— — — and — — — (Al) 
The integral of a function F(ri) then becomes 

/oo />oo />OC 1 />oo 

drF{r,) = — d V F( V ) and / dfF{rj) = — / d V F( V ). (A2) 
-oo 7" J-oo J-oo 7" ./-co 

For a double integral of the kind that appears in the formulas for the linear momentum and energy, Eq. i|14[l . and the 
angular momentum, Eq. (|17f) . we may pick r = f = 0. Then 

oo />0 i poo pO j poo p—4> 

1 ' dr) [ d4>F{r,) = -^— ( dn{-tj)F(n), (A3) 



77^ 2 V-oo Jo 77^ 2 
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and 



dr \ d?F(r]) = - 

oo JO 



1 f u Z" 00 I 
77^ L J x •/,, 



7 _ Q2 , d0 / 



77^ y Q / d <^) 



1 



77^ 2 7 



dry r)F(rj); 



adding the last two equalities, we have 



oo />0 />0 roc 

dr df- dr df ) F(rj) 

ii ./-oo J-cc JO 



77f2 2 J _ 



dr/ riF(r)). 



(A4) 



(A5) 



Our next task is to evaluate integrals of the form / drX a and / drX a , where X a and X a are the vectors 



associated with m and m, respectively. In terms of components along an orthonormal frame, the vectors have the 
form 



X a = X\rj) t a + X™{r]) m a + X*(r]) $ a + X z {r]) z c 



(A6) 



X a = t a + X™(r]) w a + X^r?) Q + X z (r]) z° 

The corresponding expressions for the derivatives are 



X a 



dX a 
~d7~ 



-70 + ^ 

\ dr? dr? 



dA"^ 



dX z 



d?7 dry 



(A7) 



(A8) 



dX z 



z a - Xf w a + X™ ^ 



(A9) 



dr \ dr] dr] dr] dr] 

where Eq. (|A1I) is used. Then the integral of I" with respect to f , with fixed r (and hence fixed </>), is then given by 

_dX a 7 



dfX a 



df- 



dr 



7 



dr) 



dx \« + d -^^ + d -^r + d -^ z « + x^«-x-r 



dr\ dr] 



dr] dr] 



=const 



r=const 



X^dr] VD a + X™ dr] x " 



1 

7 

Similarly, the integral of X a with respect to r, with fixed f (and hence fixed <fi) is 
dX a 7 



(A10) 



drX a = / dr 



df 



7 J- 



dr) 



a dX™ dX^ , dX Z v-Z-a, vvj To 

t -\ W H ; H ; Z —X V ZD +X <p 



dr] dr] dr] dr] 

- [1°]^+ / dr] (-X^ m a - X™ 4> a 

X*dT] w a + { X™ dr] 



r=const 



f=const 



(All) 



Note that basis vectors w a and 6 a are functions of d>\ and n7 Q and 6 a are functions of c 



23 



2. Integral formulas for half-retarded + half-advanced Green function 



Integrals involving the half-retarded + half-advanced Green function S(w) are derived in this section. As mentioned 
in Sec. IIII CI the solutions to w(rf) — are r\ — it ± ip. It is assumed that any contribution from r/ = ±oo vanishes. 
From Eq. we have 



dw\ 1 



n 2 



2 ip + vv sin ip ' 



and hence 



/oc 
dr,S(w)F(r,)= ]T 
-OO ^_^_L, 



dw 



rj—TT^ip 

For integrals involving a derivative of S(w) 



dr) 



n 2 



i 



2 (^9 -f- vv sin 93 



{^(tt + ip) + F{tt - <p)} 



d 

di] 



dw 
dr/ 



F^) 



Using Eq. (|A12|I and its derivative, 



d ( dw \ 1 

dn \dn) v=7r±if 



the integral becomes 



dn5'iw)F(n) = 



n 4 



1 



il 2 1 + vv cos ip 
2 ((p + vv simp) 2 

I + vv cos ip 



dF. , dF, 
4 + iro sin ip) 2 i dn dn ip + vvsmip 

Finally, the following relation is useful for computation of Eq. (|78|) in Sec. IIII DI 



[F(tt + <p) + F{n -ip)]\. (A16) 



dn s'(w) wF(n) 



dnS(w) Fin) 



dn S(w) w 



where, in the second term of the first equality, a derivative 
larger than 0(w _1 ) ) in the neighborhood of w = 0. 



dn 



dn 

dw 



dn 
Fin) 



dn 

dw 



Fin) 



dn^iw) Fin), 



(A12) 



(A13) 



(A14) 



(A15) 



(A17) 



is assumed to be nonsingular (has an order 



APPENDIX B: POST-NEWTONIAN CORRECTIONS TO E AND L 

In presenting relations used to compute corrections from Apn to E, L, and the equations of motion, we avoid 
repetition by omitting barred <-*■ unbarred versions. The corrections to the parametrization-invariant post-Minkowski 
action are calculated, assuming the parametrization 152(1 . We begin with the simpler correction term, (|124fl . Using 
the relations from Sect. IIII O with R in the form -R 2 = yp- (v 2 + v 2 — 2vvcosn), we obtain, 

mfh(m + fh) , „. 1 , . . ns 1 
A PN = -5(t-t) y — '-i-x a x a )^i-x x^ 

. r2 3 mfhim + fh) 

= — (77 — 7Tj - — = . (Bl) 

2 v 2 + v 2 — 2vv cos n 

Derivatives appearing in the equations of motion and the conserved quantities have the form 

<9A PN , £1 4 mmim + rh) 

- S fa-Tr)— - , a2 — , (B2) 



d it — t) ' 2 v 2 + v 2 — 2vv cos n ' 

9Apn 4 mfh(m + fh) 



s„4 mm(m + m) f. _ . _ -r 1 

g(7? -tt)S2 — \iv-v cos n)^a - v smnfia \ , (B3) 

a(x a — a; a ) (ir + v 2 — 2iwcos7/)- ;: I J 

— = %-7r)— 2—^2 ^ ) (B4) 

ox p 2 v 2 + v 2 — 2 w cos n 

dApN _ dA PN _ _ A} 3 mfhjm + fh)j 

dx° ~ dx a ~ [V n) 2 v 2 +v 2 -2vvcosn' [ ' 
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The index a of x a is spatial. The post-Newtonian correction term to the tu-component of the equation of motion is 
then 



df 



<9ApN d <9Apn 



dR a dr dx a 
The correction to the angular momentum is 

1 
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epN — SI^pn — — 



^ftc" 7" 



77$7 2 



mm(m + to)£1 2 



77 



j —i 



ri?7 < f5(ry — n 
{ 

d6(r) — 7r) 1 



ww 77 sin 77 



w 2 + w 2 - 2ww cos 77 (v 2 +v 2 -2vv cos riY 



+ - 



d?7 2 w 2 + w 2 — 2ww cos 



i} 



mm(m + m)Q 
277 (w + v) 2 



1 



-O ^pn- 



Finally, 



(B6) 
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e PN = -Q4>N = V ; 2 . B7 

2 277(w + w) z 

The analogous corrections from the special-relativistically covariant post-Newtonian correction are as follows. 
Correction to the equation of motion: 



J —1 
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where a = ±1 for 77 — tt — ww sin 77 ^ 0, respectively. 
Correction to the angular momentum: 
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Correction to the energy: 
6spn — f^SPN = 
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mm(m + to) 
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Finally, 
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APPENDIX C: SCALAR INTERACTION 

A scalar field is described by a Fokker action of the form (JSJ , with 



A = qqS{w) (~i a x a ) 1/2 (-5 Q i tt ) 1/2 



(CI) 



The equation of motion for two particles in circular orbit has, for to, the form 

d mx, 
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with the corresponding barred<->unbarred equation for to. 3 

The formalism developed earlier quickly leads to expressions for angular momentum and energy: 
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mx a k a + — 
U7 



, dx a 



mx a k a + — — 
SZ7 



2 9 

— n 2 

77 



dfjr]^R a k 
ow 



m7 (l — u 2 ) — -^7 / dr}5(w) r y (l — u 2 ) + TO7 (l — w 2 ) — — — / drj8{w)^i (l — u 2 ) 
"7 J-00 "7 



77$7 3 y_ c 

m to 2gg 
7 7 77$7 



dryi5' (10)77 (r/ — n — vv sin rj) 



drjS(w) 



77SI 



dr]5'(w)rj— — 
drj 



3 An equivalent solution has been obtained independently by Jean-Philippe Bruneton and Gilles Esposito-Farese l23l . 
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Then 
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APPENDIX D: DERIVATIONS OF THE MOMENTUM AND ANGULAR MOMENTUM FORMULAS 

In this appendix, we present calculations and useful expressions relating to the momentum and the angular mo- 
mentum formulas, for an arbitrary interaction term A(x — x, x, x). 



1. Derivations of Eqs. (1141) and (1171) 

Dettman and Schild p"s| derived the momentum and the angular momentum formulas from Poincare invariance 
of the general Fokker action. We follow their calculation for the case of our action for two point particles to derive 
Eqs. (fHfl and l(T7|) . from the varied action ©. 

We first rewrite Eq. (JSJ) for 51, substituting the equation of motion for each particle, Eqs. (|10fl and 
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(Dl) 
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dr dr 

We next restrict the variations of the trajectory of each particle, 5x a and dx a , to infinitesimal Poincare transfor- 
mations (which we will subsequently take to have the forms (112t and I^Jl). Since the Fokker action, and hence the 
interaction A, is Poincare invariant, the interaction term A satisfies the identity 



A(ir. + 8x — x — 5x, x + Sx, x + 6x) — A(x — x, x, x). 



(D3) 
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Expanding the identity to first order in Sx and Sx, we have 



<9i? Q <9i? Q <9± Q cte" 



(D4) 



Using this relation, the variation of the action I(D1(I becomes 
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The double integral in Eq. i|D5(l is rearranged to separate the contribution from (tx,t\) and (r2,T2) explicitly as 
follows, 



Tl •/ — OO 



— OO \>/Tl 



ri Jr 2 
00 /"fa 



Tl W — OO •/ — OO 
T2 /"Tl 



00/ ^ r 2 

OO />Tl />Tl 



(D6) 
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Finally, the variation <57 with respect to an infinitesimal Poincare transformation takes the form 



SI(t 1 ,t 2 ,t 1 ,t 2 ) 



mx 



{—x 1 x~ 1 ) 2 

OC t>T 



00 , dA 
dr- 

—00 





T2 




+ 







5x Q 



CX) 



dR a dx a 



drdr 
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To derive the conserved momentum, Eq. 114(1 , a translation by a constant vector 1(12(1 is substituted in l(D7(l : 



SI 

Sa a 



(-i 7 a;T)2 



_9A 



OO />T 



— OO J —OQ J T 



dA 
OR"' 



, <9A 
dr—r- 

(-x 7 x-y)i J-00 dx a 



drdr 
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where we used Sx a = da a / dr = 0. 

The angular momentum similarly corresponds to an infinitesimal rotation 1 (15(1 : 



m (x a if3 — xpx a ) 
(-XjX'y)i 

ffl (x a X(3 — XfiXa) 
00 />r 



— 00 J —00 J r 



dr 1 x a 
V 


dA 


di? 


dr 1 x a 


dA 


d¥ 


<9A 




Xa dBP 


- Xf3 





dA \ 




~ Xf3 


di^J 








Tl 






f% 




dA \~ 




~ Xf3 


dx a ) 








Tl 


dA 


dA 


dR a 


+ Xa d±p 



X0 



dA 

dlb* 



drdr 



(D9) 



where Sx a — e a ^ifj is used. 

An analogous derivation of the momentum and angular momentum for the affinely parametrized action results in 
the replacements 



(— x 1 x'<) i - 



( — X~yXl): 



(D10) 
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in the momentum equation i|D8|l . and 

m (x 



( — XyX"/)! 



in the angular momentum equation i|D9|) . 

The double integrals in Eas. (|D7|) - (|D9(l have a useful alternative form: 



I .x m(x a X0 XpXa) -i- - - - \ fT\-\-\\ 

m(x a xp - X(jX a ) , 7 7 ,i > m(x a x l 3 - Xf3X a ) , (Dll) 



drdr = 



dr 



drdr — 
fir - I df 



drdr 



dr. 
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The terms depending on the proper time of each path r and r can be separated by writing the double integral in the 
above form, i|D12(l . as seen in Eqs. I|D18|) and i|D20|) . These are used in the proof of the first law in Sec. IIVI 



2. Formulas for E, L, and E — Q,L 

Next, we derive expressions for the nonzero components of the angular momentum L := £12(1", f), the 4-momentum 
E := —P a (r,f)t a , and a combination E — VtL for parametrization invariant action. The formulas for the affinely 
parametrized action can be obtained by replacements analogous to Eqs. 1D10|) and 1D11|) . which are not repeated. 

We begin with the angular momentum: Since the basis </>" and <p a at positions of particles {m, x} and {m, x} have 
components in Cartesian coordinate (f> a — (— a; 2 , 2; 1 , 0) = (— x 2 ,xi,Q) and 4> a = (—x 2 ^ 1 ,®) = (— x 2l x~i, 0), one has 
relations 



X1A2 — X2A1 = A a (f> a and X1A2 — x 2 A\ = A a <fi a , 



(D13) 



where A a and A a are vectors at positions m and m. Applying the relations l|D13|) to the Li 2 component of the 
angular momentum equation i|17|) yields 



L = 



(-X/3X 



00 , <9A 
dr—— 1 

dx a 



— OC J T 



(r) 
dA , 



dR a 



{-XfjXl 3 ) 

dA 



T + 



drdr, 



dr 



dA 

dx a< 
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where we define <j) a 
corresponding form 



-£f<fi a - From the definition of the linear momentum Q14[l. the energy E = —P a t a has the 



mx a t a 



E 



The combination E — flL then becomes 



, dx a 



dA 



dR a 



(r) 



t a drdr. 



mx a t a 
{—x~ l x 1 Y- 



dr-^t° 

dx a 



(r) 
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E-VLL = - 



dA« 

, dx a 



(r) ~ 



mx a k a 



dA° 
dx a 



( + g-np ) drdr. 
\dR a dx a Y 1 
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One can write the energy l|D15(l and the angular momentum l)D14|l in a form related to the one particle energy and 
angular momentum defined in Sec. IIVI Eqs. 1)11)4)) . and l|lU(i|) (|1U8H . together with the relation for the double integral 
(|D12p . Interestingly, the resulting formulas for the energy and angular momentum depend separately on the proper 
time of each path, r and f. Using the property 



dA _ dA 
dR a ~ dx" 



dA 
dx 7 ^' 



(D17) 
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we can rewrite the energy (|D1 51 in terms of the one-particle potential U m and Um, 

E = 



mx a t a ^ d 
(-i T iT)5 dx a 



poo 




/ dfA 


(r)- 


1 —oo 





mx a t a t<x d 



drk 



(?) 



dr 



df 



mx a t a dU m 
k dx<* 

g dlifh, 

dx a 



(—XjX^) 

fhx a t a 



df 

O 

(r)- 



r/r 



<9A 



dri"^ 



dft° 



. dU rn 

dx a 
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To relate the angular momentum (|D14J| to one-particle angular momentum, we notice that the relation (|D4J| . 
together with properties HD13|) and (|D17|) . implies 



; 5A 
dx a 



dk 



+ < 



, dk 
dx a 



dk 



0. 



dx a dx a dx a 
Like the energy, the angular momentum can be written in terms of the potentials U m and U rTl 



(D19) 



L = 



(-xpx 13 ) 
dr 



dx a 
df - 



h +<t> 9^ 



dfA 
df 



(r) + 

O 

dr 



T+' 



(-i i ai' 3 ) 

fhx a (j) a 
(-cepx 13 ) 



T + 



- dUm 

dx a 



(r) 
(?) 



(-i/ji^) 



d 
da; 01 



drk 



(?) 



df 



dx a 
dx a 



dx a 

■- dUfa 

dx a 
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The terms in the brackets in the above Eqs. 1|D18|) and (|D20(I are the one-particle energy and angular momentum; 
the contribution from each particle which is moving in the field produced by the other particle (see Sec. IIV|1 . 

It is also noticeable that the formulas for E, L and E — flL, Eqs. (|D15|I . (|D14|I . and (|D16|) . can be written in a 
common form, because of a property t a := dt a /dr — 0, k a = Q<p a accordingly, and those for corresponding barred 
quantities. Writing Q to represent these conserved quantities, E, L, and E — ilL, and C, a for the associated vectors, 
i Q , #\ and k a , we have 



Q = 



mx a ( a 
{—xp±P)i 



, dx a 



mx a ( a 
(-xpxP)i 



— OO J T 



(r) 

dR a dx a 1 



, dx a 



(?) 



Also for the expressions in terms of the one-particle potentials, we have 

mx a C, a 



Q = 



{-±p±P)2 

fflX a C a 



T+C 



(-xpxP] 



T+C 



a dU m 

dx a 

dUfa 

dx a 



(r) 
(?) 



dr ( 

1 dx a 



+ C 



L a dU rn 

dx a 



dx a 



dx° 



In Eqs. (|112|) - (|115|1 in Sec. II VI the integrand of Eq. (|D22|) is written in a short form defined by 



^ dx<* " 



a dU m ^ a ^j y _^ a dU m i --: a dUm 



dx a 



dx a 



+ C 



dx a 
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As seen in the above expression (ID22() . contributions to the total energy and angular momentum from time r and 
t are separated. It means that each piece conserves independently; E L and E — QL, are a sum of two conserved 
quantities associated with r and f. 



For circular motion, the velocities are given by x a — jk a and x a = "fk a , and the parametrization-invariant 
interaction term satisfies Eq. (|99|l (Eulcr's relation for a homogeneous function of degree one in x a and x a ); we then 
have 



E-m = - 



771 X q. k 



l 



{-x p xP)^ 7 



poo 




/ dfA 


(r) - 


t — oo 





— OO J T 



dA 



k a 



mx a k a 
(-xpxf 3 )i 
dA 



drA 



tt<j) a drdf. 



dR a dx 

Only the radial components of the equations of motion, (|10|l and are nontrivial, and they take the form 



(D24) 



d_ 

dl 

d_ 

df 
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(-x 7 a;T)5 



mx 



( — XjX^) 



dA 

(XT— VJ 

, dR a 



, 5A _ c 
dr— =— rar 



7 » 2 

v 

T; 



df 



dr 



dA 

dx a< 

dA - 
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In the above we used Eqs. (|A10(1 and ijAllfl in Appendix I A II to eliminate the r derivative in the f integral (and 
respectively f derivative in the r integral), then transformed the integration variable rj back to r (f) using (|A2|) . 

The radial component of the accelerations 'x a w a = —j 2 vfl and x a -uj a = — 7 2 vSl are related to x a (j) a = jv 2 /Q and 
x a 4> a = jv 2 /Q, respectively, by 



rax 



7 f7 2 dr 



(— i 7 i 7 )^ 



07 Q and 



7 o 2 rff 



17M e 



( — X 1 X~ I ) 2 
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where x a x a = = used. With this relation l|D27|l . the radial equations of motion l|D25|) and (|D2f3|) are 

substituted to further simplify the angular momentum formula l|D14|) : 



L 



v 



dr— tu 

8R a 



(r) 



7 ft 2 



dA 
dR^' 



00 , dA 
dr - - w 

— oo 

dA ■ 



dR a 
drdf. 
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3. Conserved quantities in Post-Minkowski gravity in terms of h a p and h at 



As shown in Sec. IIII Al the parametrization invariant interaction (12211 leads to equations of motion, Eqs. (|28|l and 
()29p. that describe point particles in a post-Minkowski approximation with half-advanced + half-retarded fields. 

The momentum and angular momentum formulas, in the form of a sum of one-particle contributions, can be written 
in terms of h a p and /i a /3- We present these formulas valid for arbitrary particle trajectories m and m. 

We begin by rewriting the variation of / with respect to an infinitesimal Poincare transformation, (|D5() . in terms 
of one-particle potentials U m and Ufa, 



8I{t 1 ,t 2 ,ti,t 2 ) 
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The conserved momentum and angular momentum are then derived as follows, 



P a (r,f) 



(— i 7 i 7 ) 

fhX n 



T + 
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X0 
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(— x~ f x' 1 )i 
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X/3 



~dx° 
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dU m 
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For the interaction of Eq. 122[l and corresponding fields l|28|) and (|29[1 . the one-particle potentials are related to the 
fields h a /3 and h a p by 



U m = 



1 



x a x 13 



clrk = -mh a fj- n 

2 {-x^y. 



Ufh, — 



1 



x a xP 



drA = —mh a R : 

2 Q VM 7 )^ 



and their derivatives with respect to position and velocity take the form 



dU m 1 dh 7 s x^x s 
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8x a 2 dx a (- Xe x*)i 
dUm 1 _ dhjs x^x 6 
~8x^ ~ 2 m 0£°"(-i 6 i«)j 



dx a 
dx a 



1 



T-a/3 



T-a/3 



h"y§ x"~^ X 
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Substituting relations (|D34|) and 1D35(1 in Eqs. (|D31(I and l|D32|) . we obtain explicit expressions for the momentum 
and angular momentum in terms of h a /s and h a /3, 



P a (r,f) 
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and 
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For completeness, we give the corresponding form of the conserved quantity Q associated with a Killing vector of 
Minkowski space, rewriting Eq. I|D22|) . in terms of h a p and h a p- Applying a property of £ Q , x a C, a = x a x"'V a Ci3 = 0, 
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we have 
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1 fl/ryfiXs^ 3(j 
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For the affinely parametrized interaction (|36|) , we obtain analogous expressions for momentum and angular momen- 
tum in terms of h a p and h a p. Since the form l|36|l relates to the solutions 

h a p(x) — 4to J d? 5{w) (^x a xp — -r] a px 1 x' y ^J , (D39) 

h a fi(x) = 4m J <1t5(w) (x a ±p - ^rjapx^x^ , (D40) 
the one-particle potentials for the affinely parametrized interaction are written, 
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drA — —m h a fix a x^ , 



Un 



1 



drA = —m h a px a x 13 , 



and their derivatives with respect to the position and to the velocity become 



(D41) 
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Replacements (|D10|I . IjDlip and substitution of (|D42|) and (|D43|) in Eqs. 1D31|) and l|D32|> yield formulas for the 
momentum and angular momentum the expressions, 
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and 



L a p(r, t) = to {x a xp - xpx a ) + (x a hpi - xph ai ) i 7 (r) 



to / dr 

' — oo 



- %#a) + {xjipy - xph ai ) X 1 (f ) 
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Finally, the charge Q has for the affinely parametrized interaction the form 



+ m(r) af 3 + hap) ( a x"(T) 



rn 
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